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Results from the ICARUS T600 module 

A measurement of the ji decay spectrum 
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Abstract. We have studied the /x decay energy spectrum from a sample of stopping /x events acquired 
during the test run of the ICARUS T600 prototype. This detector allows the spatial reconstruction of the 
events with fine granularity, hence the precise measurement of the ^ range and dE/dx with high sampling 
rate. This information is used to compute the correction factors needed for the calorimetric reconstruction. 
The Michel p parameter is then measured by comparison of the experimental and Monte Carlo simulated 
^ decay spectra, obtaining p = 0.72 ± 0.06 (stat.) ± 0.08 (syst.). 

PACS. 13.15.+g Neutrino interactions - 13.35.Bv Decays of muons 
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1 Introduction 

The sample of events in which a muon enters the detec- 
tor, stops and eventually decays in the detector's sensi- 
tive volume -hereafter called stopping muon sample- con- 
stitutes an important benchmark to evaluate the physics 
performance of the ICARUS detector. Because of their 
simple geometry, stopping muon events are relatively easy 
to reconstruct in space, which allows the computation of 
the different correction factors needed for the calorimet- 
ric reconstruction. Thus, we can study the muon decay 
spectrum and measure the Michel p parameter, which 
constitutes the first physics measurement performed with 
the ICARUS novel detection technology, and proves that 
the technique is mature enough to produce competitive 
physics results. 

Muon decay was first described in a model-independent 
way by Michel , using the most general local, derivative- 
free, lepton-number conserving, four fermion interaction. 
For unpolarized muons, the decay probability is given by: 

^(x; p, ^) = Ul -x) + ^p(4x - 3)+ (1) 



3r? 



l-x 



X = is the reduced energy (ranging from mf,/Emax to 
\)] Ee and We are respectively the total energy and mass 
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of the electron produced in the decay; Emax = 52.8 MeV 
is the end-point of the spectrum; f{x) is the term ac- 
counting for the first order radiative corrections assum- 
ing a local V-A interaction Q^; finally, p and rj are the 
so-called Michel parameters, defined in terms of bilinear 
combinations of the coupling constants of the general four 
fermion interaction, and hence depending on the type of 
interaction governing the decay process. For the Standard 
Model (SM) V-A interaction, the parameters take the val- 
ues psM = 0.75 and tjsm = 0. 



2 The ICARUS T600 detector 

ICARUS T600 121 is a large cryostat divided in two iden- 
tical, adjacent half-modules of internal dimensions 3.6 x 
3.9 X 19.9 m'^, each containing more than 300 t of liquid 
argon (LAr). Each half- module houses an internal detec- 
tor composed by two Time Projection Chambers (TPC), 
the field shaping system, monitors, probes, PMT's, and is 
externally surrounded by a set of thermal insulations lay- 
ers. Each TPC is formed by three parallel planes of wires, 

3 mm apart, oriented at 0, ±60° angles, of 3 mm pitch par- 
allel wires, positioned onto the longest walls of the half- 
module. The cathode plane is parallel and equidistant to 
the wire planes at each TPC. A high voltage system pro- 
duces a uniform electric field, perpendicular to the wire 
planes, forcing the drift of the ionization electrons (the 
maximum drift path is 1.5 m). 

The ionization electrons produced in the LAr active 
volume drift perpendicularly to the wire planes due to the 
applied electric field, inducing a signal {hit) on the wires 
near which they are drifting while approaching the dif- 
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Fig. 1. Run 966 Event 8 Right chamber: muon decay event 
views corresponding to the Collection (left) and second Induc- 
tion (right) wire planes. 



ferent wire planes. By appropriate biasing, the first set of 
planes can be made non-destructive (Induction planes), so 
that the charge is finaUy coUected in the last plane {Col- 
lection plane) . Each wire plane provides a two-dimensional 
projection (view) of the event, where the position in one 
coordinate is constrained by the hit wire, while the signal 
timing with respect to the trigger gives the position along 
the drift direction. 



3 Data selection and reconstruction 

Data were acquired during the T600 technical run 3 . The 
off-line data selection was carried out by visual scanning 
using topological criteria. A total of 5830 triggers were 
scanned, containing 4548 stopping muon events, out of 
which 3370 were further selected by the preliminary qual- 
ity cuts. 

The spatial reconstruction of the muon tracks is needed 
in order to compute the corrections entering the calorimet- 
ric reconstruction of the events. A detailed description of 
the spatial reconstruction tools has been reported else- 
where 

The ionization charge is precisely measured in the Col- 
lection wire plane. The energy associated to a given hit is 
related to the collected charge by means of 



E 



CW 
R 



= (t-*0)/Te 



(2) 



23.6tg j eV 



where C is the calibration factor 0; W 
is the average energy needed for the creation of an elec- 
tron-ion pair 0; R the electron- ion recombination factor; 
(t — to) the time of drift of the electrons; the drift elec- 
tron lifetime, which parametrizes the attachment of drift 
electrons to impurities in LAr; and Q the measured charge. 
R, to and Te are extracted from the reconstructed muon 
tracks 0|7], essentially by tuning them so that the mea- 
sured energy corresponds to the theoretical expectation 




1000 990 



960 950 
z (cm) 



Fig. 2. Run 966 Event 8 Right chamber: fully reconstructed 
muon decay event. 



for stopping muons. This method determines the electron 
energy in a bias-free way, since all the correction param- 
eters are tuned using exclusively muon tracks. For Te we 
have obtained values ranging from 1.2 to 1.7 ms (depend- 
ing on the data taking period), meaning free paths for 
drift electrons between 1.9 and 2.6 m, to be compared 
with the maximal drift distance 1.5 m. The measured re- 
combination factor is i?mip = 0.640 ± 0.013 for minimum 
ionizing particles and R~^ linearly increasing with slope 
0.11 ±0.01 cm/MeV. 

Figure n shows the initial two-dimensional projections 
of a typical stopping muon event, produced by the Col- 
lection and the second Induction planes, respectively. On 
Figure 13 the fully reconstructed event is shown. 



4 Results 

A total of 1858 electrons are available for the determi- 
nation of the Michel p parameter. The measured spec- 
trum corresponds to the fraction of the electron energy 
lost by ionization in LAr since no attempt to recover the 
energy loss by bremsstrahlung radiation has been carried 
out. We compare the spectrum with the one obtained for 
a Monte Carlo (MC) simulated event sample, containing 
10000 electron events from muon decays, generated using 
FLUKA 8^. The simulation includes all detector effects 
except for the presence of impurities, whose effects are 
evaluated from data and included in average in the final 
MC distribution [7]. The measured and simulated energy 
spectra are compared in Figure |3| They are found to be 
in good agreement {x^/ndf = 14.0/20). 

We measure p while constraining rj within its experi- 
mentally allowed interval (-0.020 < tj < 0.006) . The 
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Fig. 3. Muon decay spectra in the ICARUS T600 detector. 
The plot shows the measured distribution (dots) compared 
both to the SM expectation (solid line) and the best fit with p 
and f) as free parameters (dashed line). 

spectrum for an arbitrary pair of values, p and 77, is built 
by weighting the MC events with a factor -5^(^"C'P''') 

where xmg = E-Mc/Emax, and E'mc is the generated en- 
ergy. We extract the value of p as that for which the best 
fit between the simulated and measured energy spectra is 
obtained, which yields p — 0.72 ± 0.06, where the error 
is of statistical origin and includes the correlation with 77. 
Figure 121 shows the spectrum for the best fit of p. 

There are two types of source of systematic uncertain- 
ties that can affect this measurement, namely: an underes- 
timation of the energy resolution and a systematic shift in 
the global energy scale. With the help of the MC sample, 
the contributions to the total systematic error are esti- 
mated in ±0.01 and ±0.08, respectively. Therefore, the 
total error is dominated at this level by the systematics 
due to the uncertainty on the energy scale, which stresses 
the importance of a fine calibration of the detector. The 
final result is: 

p = 0.72 ± 0.06 (stat.) ± 0.08 (syst.) (3) 

compatible within the error bounds with the V-A value. 

The Michel parameter p from p decays has been mea- 
sured in the the late 60's by several authors ^U] with a 
precision of about 0.4%. Such results have been obtained 
using dedicated experiments involving data samples of 
typically several hundred thousand events. More recently, 
the data from electron-positron colliders have been used 
to measure the Michel parameters of the purely leptonic 
T decay These measurements have an accuracy of 
about 5% and are based on the analysis of samples typ- 
ically amounting up to several ten thousand events. Our 
new result is not competitive with those obtained from 



p decay, and barely with those obtained from t decay. 
However, it must be remarked that this result has been 
obtained using 1858 muon decay events with a non opti- 
mized experiment. This result stresses the capabilities of 
the ICARUS technology to produce robust physics results. 

5 Conclusions 

We have performed the first physics measurement with 
the ICARUS LAr TPC detection technique, the Michel p 
parameter, from the detailed study of muon decay spec- 
trum using the stopping muon event sample from the 
ICARUS T600 detector test run. We have obtained p = 
0.72±0.06 (stat.)±0.08 (syst.), in agreement with the SM 
value. This measurement involves the exploitation of both 
the spatial and calorimetric reconstruction capabilities of 
the detector. Therefore, the obtained result constitutes a 
proof of the matureness of the detection technique to pro- 
duce high quality physics results, in particular for neutrino 
physics when it will be installed in the Gran Sasso under- 
ground laboratory. 
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